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Abstract. Motivated by the Double Gyroid nanowire network we 
develop methods to detect Dirac points and classify level crossings, 
aka. singularities in the spectrum of a family of Hamiltonians. 

The approach we use is singularity theory. Using this language, 
we obtain a characterization of Dirac points and also show that the 
branching behavior of the level crossings is given by an unfolding 
of A n type singularities. Which type of singularity occurs can be 
read off a characteristic region inside the miniversal unfolding of 
an At singularity. 

We then apply these methods in the setting of families of graph 
Hamiltonians, such as those for wire networks. In the particular 
case of the Double Gyroid we analytically classify its singularities 
and show that it has Dirac points. This indicates that nanowire 
systems of this type should have very special physical properties. 
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Introduction 

In many physical situations one is led to a family of finite dimen- 
sional Hamiltonians defined over some parameter space (base) B and 
would like to consider and classify the level crossings that appear when 
one changes the parameters. We came upon this general context when 
studying the C*-geometry of wire networks in general and the Dou- 
ble Gyroid network in particular [lj|2). These networks are spatially 
periodic, and the base B is a torus spanned by the values of the quasi- 
momenta. The dependence of energy eigenvalues on quasimomenta 
determines the band structure; the simplest type of band crossing, a 
conical intersection, is often referred to as a Dirac point. Band cross- 
ings are interesting because they may be responsible for new physical 
phenomena (as well known, for instance, in the case of Dirac points 
in graphene). Of particular interest are Dirac points in triply periodic 
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materials, such as the Gyroid network: they can be viewed as mag- 
netic monopoles in the 3-dimensional parameter space [3] and as such 
are expected to be stable under small deformations of the Hamilton- 
ian. In order to study various types of level crossings, we first widen 
the context to that of a family of Hamiltonians over an arbitrary base, 
and then apply the general results to our initial problem, in which the 
base is actually the compact n-torus. 

To be more specific, we consider a differentiable map from an n— 
dimensional manifold B to the set of k x k Hermitian matrices for 
a fixed k. The case of interest will be B = T n , the n-dimensional 
torus, and it can be thought of as the space of momenta. Our results 
about the multiplicities in the spectrum are obtained using singularity 
theory [2]. They are twofold. First we give an analytic way of finding 
all Dirac points, by considering the energy levels as the zero set of a 
smooth function P on B x K. Using the Morse Lemma (Theorem 1.1) 
we explain that a Dirac point in this context and language is an isolated 

Ai singularity with the signature ( h) (or (H 1 — ) depending 

on the sign of P) for the function P. This effectively uses an ambient 
space to embed the conical singularity. 

Then by considering the energy levels as a singular fibration over 
the base space of momenta, we classify the possible singularities in 
the fibers. The fibration in question is the first projection BxR-> 
R. For this we again use singularity theory, more precisely that of 
the singularity A^_i and its miniversal unfolding. In particular, we 
define a characteristic map of the base B of the family to the base 
A = C of the miniversal unfolding. From its image, which we call 
the characteristic region, one can read off many details, such as at what 
points degeneracies occur, and what their nature is. Degeneracies occur 
precisely over the points of intersection of the characteristic region with 
the discriminant locus. 

This point of view allows us to classify the possible degeneracies as 
those appearing in the discriminant locus or swallowtail of the A^-\ 
singularity. These are known by a theorem of Grothendieck [5] on the 
singularities appearing in the fibers of the miniversal unfolding associ- 
ated to the singularities corresponding to Dynkin diagrams. Namely, 
these are precisely those obtained by deleting vertices (and all incident 
edges) of that of A^-i and hence are of the form (A ni , . . . , A nr ) for suit- 
able rii. This corresponds to simultaneous crossing of n\ + 1, . . . , n r + 1 
levels. How these levels cross or equivalently how the singularity un- 
folds is encoded in the characteristic map and can qualitatively be read 
off from the characteristic region. 
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One way to view this result is as a more precise and more general 
version of the von Neumann- Wigner theorem [fj]. Indeed for the full 
family of traceless 2x2 Hamiltonians in their standard parameteriza- 
tion over M. 3 , we reproduce that the locus of degeneracy is of codimen- 
sion 3. More precisely, there is only one point G M 3 in the preimage 
of the characteristic map restricted to the discriminant. 

If the family is more complicated however, our methods tell us where 
degeneracies can occur and how the levels cross. In this case the codi- 
mension 3 is not universally true any more. Among the graphs we 
study, we exhibit families, where the codimension of the degenerate 
locus is 2, 3 or 1. The precise dimension count comes from the in- 
tersection of the discriminant with the characteristic region and the 
dimension of the fibers of the characteristic map. For isolated singu- 
larities such as Dirac points one needs that the particular fiber of the 
characteristic map over the corresponding point in the discriminant is 
O-dimensional. This is true in the von Neumann-Wigner case and this 
fact corresponds to the "extra equations" as we explain. 

Our main aim of application is the commutative and non-commutative 
C7*-geometry of wire networks [TJ[2] in their description as graph Hamil- 
tonians. The input data for the general theory are a graph T, which 
is embedded in R n , the crystal graph, together with a maximal sym- 
metry group L isomorphic to Z n and a constant magnetic field 2-form 
given by a skew symmetric matrix O. A fundamental role in the whole 
theory is played by the abstract quotient graph f := T/L. This graph, 
together with the induced data of the magnetic field and the embed- 
ding, was used to define the Harper Hamiltonian and the relevant C* 
algebra SB. This algebra which we called the Bellissard-Harper algebra 
is the minimal C*-algebra generated by the magnetic translations cor- 
responding to L and the Harper Hamiltonian H. In [2] we showed that 
SB embeds into M/^Tq), the k x k matrices over the non-commutative 
torus Tq. Here O contains the information of the 5-field for the lattice 
T and k is the number of vertices of T which is the number of sites in 
a primitive cell. 

One intriguing aspect is that this description is very useful even 
in the commutative case, that is in the absence of a magnetic field. 
The C* approach yields a family of finite dimensional Hamiltonians 
parameterized over a base torus T n = S 1 x • • • x S 1 . Namely, in the 
commutative case O = and Tq is the C* algebra of complex valued 
continuous functions of the n-torus T n = S 1 x • • • x S 1 . In standard 
notation Tq = C(T n ). Likewise, using the Gel'fand-Naimark theorem 
the Bellissard algebra SS is also the C* algebra of a certain compact 
Hausdorff space X; g$ = C(X), which as we showed in [l] is a branched 



4 RALPH M. KAUFMANN, SERGEI KHLEBNIKOV, AND BIRGIT KAUFMANN 



cover over T n . Physically the base T n parameterizes the momenta and 
the space X is given by the energies of H as these momenta vary. In 
this sense they give the energy bands of a one-electron system. 

The central questions that arise are the following: At what points 
do we have degenerate eigenvalues in the spectrum and which of these 
points are Dirac points? And, can these be read off from the graph T 
and its decorations, such as a spanning tree or weights with values in 
some Tq? 

Given a specific graph with C*-algebra valued weights on the edges, 
this is done by analyzing the function P and the characteristic map 
from the base torus T n to the miniversal unfolding of the A^-i singu- 
larity. We apply these considerations to the cases of the Double Gyroid 
wire network which was our initial interest and, to illustrate the con- 
cepts and the possible behaviors, we consider several other examples 
along the way including the wire networks obtained from the double 
versions of the P and D surfaces as well as the honeycomb lattice. 

Our main result here is an analytic proof that the spectrum of the 
Gyroid has four singular fibers, two of which are A2 singularities and 
two of which are of the type (Ai,Ai). We furthermore show that the 
latter two are Dirac points. This is the first analytic proof of this fact. 
The relevant family of Hamiltonians also arises in a different context [7] . 
There the authors found numerically that the singular points lie on the 
diagonal of B = T 3 (viewed as a cube with opposite faces identified) 
and obtained the spectrum on the diagonal, see also [3 A and [§]. We 



note that although this shows that on the one-dimensional sub-family 
of Hamiltonians given by the diagonal, there are two triple degeneracies 
and two two-fold double degeneracies, using this information alone one 
cannot conclude how the singular structure extends onto the full 3- 
dimensional torus. Our method gives this extension. In particular, it 
allows us to show analytically that there are no degeneracies anywhere 
else in T 3 . 

Applying our program to the honeycomb lattice, we rediscover the 
well-known Dirac points of graphene [9]. Thus it can be hoped that 
the Dirac points of the Gyroid will give rise to important new material 
properties. 

Another natural question in the case of graph Hamiltonians, which is 
addressed in a separate paper (8] , is if there are symmetries that can be 
derived from the graph setup which explain the degeneracies. The short 
answer is that (a) the symmetry group must be extended beyond the 
permutation symmetry of the graph to include phase transformations 
on the vertices ( "re-gaugings" ) , and (b) in all the wire network cases 
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corresponding to the 2 and 3 dimensional self-dual graphs: P, D, G 
and honeycomb these symmetries force all the singularities. 

1. Singularities in spectra of families of Hamiltonians 

1.1. Singularities. We briefly recall pertinent definition of singular- 
ity theory In this theory one considers germs of smooth functions 
/ : C™ — > C with critical point at with critical value up to equiva- 
lence induced by germs of diffeomorphisms. That is the germ (/, 0) is 
equivalent to a germ (/', 0) if there exists a germ (g, 0) of a diffeomor- 
phism g : C n — > C n with g(0) = such that / = fog. A singularity is 
an equivalence class of such germs. The germ (/', 0) is then also called 
the pull-back under g. 

Analogous definitions hold in the case of complex singularities, that 
is for germs functions / : C™ — > C with diffeomorphisms replaced by 
biholomorphic maps. 

A deformation of a germ / : C n — > C with base A = C fc is the 
germ at zero of a smooth map F : (C n x C k C, 0) which satisfies 
F{x, 0) = f(x). A deformation F' is equivalent to F if there is a 
smooth germ of a diffeomorphism g : C n x C k — > C n x C k at zero with 
g(x,0) = x. 

Given a deformation F with base A a smooth germ 9 : (C r , 0) — > 
(A, 0) induces a deformation F' via pull-back. F'(x, A') := F(x, 9(X}). 

A deformation F of a germ / is versal if every deformation F' of / 
is equivalent to a deformation induced from F . It is called miniversal 
if A is of minimal dimension. 

Again one can replace C with K. Also, one can pull-back not only 
germs, but actual pointed families by the same procedure. Here the 
base spaces C k are then simply replaced by smooth manifolds B. This 
yields the same local theory. 

1.2. The spectrum as a zero locus. The basic starting point of our 
analysis in this section is that given a smooth family of Hamiltonians 
over a base B the spectrum as functions on B can be alternatively 
given as the zero locus of a single function P on B x M. The map 
that associates to a point b in the base B the Hamiltonian H(b) is a 
smooth map H : B — > Herm(k), where Herm(k) are the Hermitian 
k x k matrices. Now since the matrices depend differentiably on the 
parameters, varying these eigenvalues gives rise to a cover tt : X — > B. 
Here a priori the cover is just that of sets, but one can quickly show 
that this is a cover of topological spaces, e.g. using C* geometry, see 



6 RALPH M. KAUFMANN, SERGEI KHLEBNIKOV, AND BIRGIT KAUFMANN 



The inverse image of a point B under it is the set of eigenvalues 
of H(b). As these are the zeros of the characteristic polynomial of 
H(b), we obtain another description of X as a subset in an ambient 
differentiable manifold as follows. 

Consider the trivial cover BxC-^B and its real part B x R — > B. 
On the space B x C, we consider the function P : B x C — > C given 
by P(b,z) :=det{zId-H(b)). 

The zero locus of P is exactly X. We chose the normalization, so 
that P starts with +z k . Fiberwise for the zero locus, we just get the 
eigenvalues of H(b), and since H(b) is Hermitian, we know that we have 
real eigenvalues and hence X is also the zero locus of P(z, b) contained 
inBxi. 

The description above gives X as a singular manifold. Around points 
of X at which P does not have a critical value 0, P _1 (0) is a smooth 
manifold. The points at which P is critical with critical value are 
singular. 

Physically P _1 (0) are just the energy levels. A level crossing can 
occur only at critical points of P with critical value 0. Namely, if we fix 
(6,2) 6 Bxl and is a non-critical value in a small neighborhood U of 
(6, z) then P\u\0) is a smooth manifold. More precisely, if ^ 7^ the 
implicit function theorem states that for (6 , zq) E U with P(bo, z ) = 
there exists a function z = E(b) such that P(b, E(b)) =0 and the graph 
z = E{b) is the component of the smooth manifold P|^ 1 (0) containing 
(b ,Zo)- In other words: E(t) is the dispersion relation. 

Our main results describe the singular locus of the singular manifold 
X. There are two approaches we will take. First we can look at the 
singularities of X locally where we regard X as embedded in B x R. 



We use this to find Dirac points, see S 1.3 Secondly, we can look at the 



cover 7r : B x M — > B restricted to X, that is n : X — > B. The upshot 
is that locally around a singular point x, X is the deformation of the 
singularity of the fiber over ir(x). By Grothendieck, locally in a fiber the 
only singularities that can appear are of type A n with n < k — 1. This 
point of view lets us classify these singularities and their deformations 



by means of the miniversal unfolding of the A^-i singularity, see { 1.3.1 



1.3. Dirac points as Morse or A\ singularities. From the point 
of view of material properties, one of the most interesting singularities 
that can occur are Dirac points. In the general terminology, a Dirac 
point is a conical singularity in the spectrum when two levels cross such 
that the dispersion relation is linear. 
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In our situation, this can be formalized in order to yield analyti- 
cal tools to find and classify these points without solving the Eigen- 
value equations. Instead of actually finding an expansion, we will use 
a smooth ambient space to characterize Dirac points using singularity 
theory. 

The standard cone is of the form z 2 = ^2i=i^i- Here as above we 
fix the sign of the coordinate z to be positive. We can rewrite this as 
F(t,z) = where F(t,z) = z 2 — 52»=i^i- The characteristic features 
of the function are that / has a critical point at = (0, . . . , 0) with 
critical value 0. Moreover the Hessian, i.e. the matrix of its second 

derivatives, is a quadratic form with signature (— • ■ h). In particular 

its determinant hess = det(Hess) ^ 0. In general if / has a critical 
point and hess ^ the critical point is called a Morse critical point. 

The most pertinent theorem about Morse critical points is the Morse 
Lemma. 



Theorem 1.1. [10] In a neighborhood of a nondegenerate critical point 
p of a smooth function F : M — >■ R from an n-dimensional manifold 
M there are coordinates x$ centered at p, such that in these coordinates 

F = -x\ - x\ x\ + x\ +l + h x 2 n + f(p) 

where, A is the index of the critical point. 

Now we see that a Dirac point as a germ is equivalent to the germ 
of (/, 0) above, whose index A is n — 1 or —1 if one switches the sign 
of F, that is if one regards — F = 0. Let us for the moment assume 
that the sign of F is chosen such that the signature is (— ■ ■ ■ — V) in 
the order of coordinates above. 

Such a germ is the pull-back under some diffeomorphism on the 
ambient space and the cone itself is the zero set of the function F. 
Notice that the characteristic properties of being a Morse critical point 
are invariant under the diffeomorphism. The signature has a geometric 
meaning. It says that the cone opens up on the The dispersion 

relation in the new coordinates is just the pull-back and hence also 
linear. 

In our particular case, that of the fibration 5x14 5 the function 
F = P and the role of x n should basically be that of the coordinate z 
on the fiber B x R — > B and (x±, . . . ,x n -i) should correspond to the 
variables on the n-dimensional base, so that we indeed get a physically 
sensible dispersion relation z = E(b). Since we are working with germs, 
we do not distinguish between a local neighborhood in B and the image 
of its local charts in R n . 
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Now the cone has the right orientation as long as coefficient > 0, 
as this states that the z-direction lies in the positive part of the cone. 
The dispersion relation can depend on the direction that is the cone 
might be deformed and tilted. The tilting of the cone is given by the 
partial derivatives . The cone is untilted if S-Jj- = for the base 
coordinates fej. Rephrased, we need that T(B x I) = TB © TR is a 
decomposition into a negative definite and a positive definite subspace. 

Recall that a stabilization of a singular germ /(z) is a function 
/(z, w) ± w\ ± • • • ± w 2 m . Thus the Morse singularities are then stably 
equivalent to the A\ singularity given by the germ f(z) = z 2 . 

1.3.1. Singularity Characterization of Dirac Points. Therefore, we get 
the Dirac points in the spectrum are precisely the critical points with 
critical value of P which are stabilizations of an A\ singularity with 

signature ( — h) or (H H — ) in the coordinates (b\, . . . , b n , z) 

such that TiB x M) = TB © TR is a decomposition into a negative 
definite and a positive definite subspace. We can now take this as their 
definition. 

Practically this means that we have to simultaneously solve the equa- 
tions P = 0, VP = and then check hess ^ and moreover check that 
the signature is correct, by computing the principal minors and check 
that = 0. 

1.4. The spectrum as a pull— back from the miniversal unfold- 
ing of the Ak-i singularity. 

1.4.1. The miniversal unfolding of Ak_i and the swallowtail. The 
singularity is the singularity defined by the function f(z) = z k , which 
has a critical point of order k — 1 at 0. Its miniversal deformation [I] 
is 

(1) F(a, z) = z k + a k „ 2 z k - 2 + ■ ■ ■ + a 

According to the general theory, the dimension of a miniversal defor- 
mation coincides with the dimension of the Milnor ring C[z]/(f) (also 
called the Milnor number) and the terms which are added to / are 
in 1-1 correspondence with the vector space basis (1, z, z 2 , . . . , z k ~ 2 ) of 
this ring. 

The geometry of the situation is very similar to the cover 7r consid- 



ered introduced in { 1.2 In particular, the function F(a, z) is a function 
on C k ~ l x C. Let Y := {(a,z) : F(a, z) = 0} C C fe_1 x C and consid- 
ering the trivial bundle C fc ^ x C 4 C fc_1 . A gain we get a branched 
cover 7r : Y — > C fc_1 . The inverse image under ir is the set of roots 
of the polynomial. Generically there are k of these roots. However, 
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Figure 1. Zero locus in the A 2 and A 3 singularities 



over a subset D C C of the base space the number of inverse images 
drops as there are multiple roots. This set is known as the discrim- 
inant locus, the swallowtail or the level bifurcation set and has been 
extensively studied (see |4j|TT]). It is the zero set of the discriminant 
of the polynomial F(z) := F(z, a) which is considered as a polynomial 
with arbitrary coefficients. The discriminant is a simple polynomial in 
the a,- and its zero set has codimension 1 
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Pictures for this locus in the A 2 and A 3 are shown in Figure [TJ The 
surface D in the A 3 singularity is known as the swallowtail. In general 
D goes by the names of discriminant, level-bifurcation set or also again 
as the swallowtail. 

1.4.2. The spectrum as a pull-back. Consider the trivial bundle B x 
C — > B as before and let P(b, z) be defined as before. We expand 



(2) P(b, z) = z k + a k ^{b)z k - 1 + a k . 2 (b)z k - 2 + ■■■ + a (b) 

where now the coefficients a« : B — > R are real since the matrix H(b) is 
Hermitian. This has great similarity to Eq. ([T|, except for the second 
leading term not vanishing. 

However a simple invertible smooth transformation s : z i— )■ z — 
a k-i/k yields a polynomial of this type. In our setup, the trans- 
formation s gives a diffeomorphism g : B x WL — > B x M. given by 
(6, z) i — y (b,z — ak-i(b)/k). This gives an equivalence between P and 
P := P o g. Now P expands as 
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(3) P(b, z) = z k + a k - 2 (b)z k ~ 2 + ■■■ + o (6) 

Let H = (oq, • • • , Qfe-a) : B — > K fe_1 , then the miniversal unfolding 
F pulls back via 5 to the deformation : B x C — > C given by 

E*(F)(b, z) = P(a(b), z). In other words P is the deformation induced 
from the miniversal deformation F via S. 

We define the characteristic map to be the coefficient map H : 5 — >■ 
C fe_1 and the characteristic region R to be the image of 5. 

Let H*(y) C 5 x C be the zero locus of then we get a pullback 

of the map 7r: S*(7r) : S*(Y) — >• £? by restricting the projection. 

Notice that if the Hamiltonians H(b) are traceless then = and 
P = P. Summing up we have the following: 

Theorem 1.2. The branched cover X B is equivalent via g to the 
pull back of the miniversal unfolding of the A^-i singularity along the 
characteristic map 5. 

Moreover if the family of Hamiltonians is traceless, the cover is the 
pull-back on the nose. 

Corollary 1.3. The possible degeneracies are pull-backs of those which 
appear in the miniversal unfolding of Ak-i- Moreover, these singulari- 
ties are present over the fibers over the real part of the discriminant. 

Since Ak-i is a simple singularity the following theorem of Grothendieck 
applies. 

Theorem 1.4. \3$ The types of singularities which appear in the swal- 
lowtail of a simple singularity are exactly those corresponding to the 
Dynkin diagrams obtained by deleting vertices and all the edges inci- 
dent to these vertices in the Dynkin diagram of the original singularity. 

If the resulting diagram is disconnected this means that there are 
several critical points with critical values in the fiber. In fact there is 
a stratification of the swallowtail S into strata . . . , X n ) according 

to the types Xj. The above theorem then states which strata are non- 
empty. 

If we take the Af. Dynkin diagram and delete I points, there are at 
most I + 1 connected components all of type A r for some r and the 
sum of their Milnor numbers r is at most k — I. Deleting points at the 
edges or next to each other, we can make the individual Minor numbers 
smaller. 

Corollary 1.5. The only possible types of singularities for nonloop 
graphs in the spectrum are (A ri , . . . , A rs ) with r« < k — s. 
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Example 1.6. In the unfolding of the A3 singularity, we have an A3 
singularity at the origin, we have A\ singularities along the smooth part 
of the swallowtail corresponding to deleting two points of A3. Along 
the cusps of the swallowtail there are A2 singularities corresponding 
to the triple degeneracy of the roots corresponding to the right two 
vertices and the left two vertices and over the double points there are 
(Ai,Ai) singularities corresponding to deleting the middle point of A3. 

Corollary 1.7. Near any singular point x G X of P there is a neigh- 
borhood of x which is a deformation of an A r singularity for some 
r < k - 1. 

Proof. If x is a singular point, then y = E(ir(x)) lies on the discrim- 
inant. Picking a neighborhood of y and pulling it back to B and to 
X, we obtain the desired deformation by restricting to the component 
that x lies in. □ 



1.5. Characteristic region. Since the families we consider are given 
by Hermitian matrices, they have positive eigenvalues. Hence, if disc 
is the discriminant function of F, then disc o S > as this function 
is is quadratic in the differences of the eigenvalues. Thus, we get that 
the characteristic region is contained in the locus of A over which the 
discriminant is non-negative. 

Notice that if B is compact connected, then the image under S of 
B is compact connected. This will be the case for graph Hamiltonians 
where B = T n , hence it will then also lie in the closure of a component 
of the real part of R n_1 \ D n i? n_1 over which the discriminant is 
positive. 

Since the discriminant is a simple polynomial, the sign changes when 
crossing the discriminant locus. This entails that the intersection of D 
with the characteristic region is only along its boundary Bd(R). Thus 
if p is an interior point of R then all the fibers of X over the inverse 
images of p under S are non-singular and moreover there is a whole 
non-singular neighborhood of each fiber. 

The fibers of X that are singular all lie over points b whose image 
sits in D n R C Bd(R). Thus the characteristic region gives 
a useful aid in studying the singularities that occur. If R fl D — 
then there are no singularities. In low dimensions this is also a great 
visualization tool. 

Notice that if B is connected, isolated intersection points only occur 
for constant maps, which means that there are no degrees of freedom. 
This cannot happen in the crystal/wire case. 
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If there are non-isolated intersection points, then their inverse images 
are singularities. Namely, by the above, there are nearby fibers, where 
the number of pre-images is k. 

Over a fiber in a component with an A r singularity, we hence know 
that there are transversal directions, where one Eigenvalue splits into 
r + 1 different pre-images. Or going into the singularity r + 1 energy 
level coalesce. Thus over a point in the stratum (A ni , . . . , A ni ) there 
are k crossings of n\ + 1, . . . , ni + 1 levels, respectively. 

1.6. Consistencies and necessary conditions. There are several 
consistencies which one might exploit for analytic or numerical solu- 
tions. Any pre-image b of a point p in the boundary of R has to satisfy 
that Js(6)j the Jacobian of H at b, does not have maximal rank. This 
is of course not a sufficient condition. We know that to get a singular 
point we need that Js(b) = 0. This again fits well with the fact that 
the discriminant restricted to R is > 0. So that disc o S has a zero 
Jacobian. 

In order for an isolated singularity at x G X, such as a Dirac point, to 
occur, a necessary condition is that the fiber H~ 1 (H(7r(a;))) is discrete. 
This takes care of the vertical direction, but of course there should 
also be no curve through x transversal to the fibers mapping to the 
discriminant. All these types of behaviors can be found in the examples 
we give. 

Another nice consistency check is given by the discriminant of P 
considered as a polynomial in z. Since any singular point x in the 
spectrum is a critical point with critical value zero, we must have 
P{x) = P z (x) = 0, which means that discriminant disc(P)(x) = 0. 
Denoting the discriminant of the A^ singularity by disc as well, we 
have disc(P)(x) = disc(E(ir(x))) = 0. 

1.7. Standard von Neumann— Wigner Example. We consider the 
family of Hamiltonians if (a, b, c) = aa x + bo y + ccr 2 , where a x , cr y , o z 
are the Pauli matrices. This gives us a family with base M 3 . It is the 
full family of traceless Hermitian 2x2 matrices. 

The usual interpretation of von Neumann- Wigner is that for a sin- 
gle level crossing one can reduce to this 2x2 family. However, this 
is basically only true in a "generic" or abstract setting and not for 
any arbitrary particular family. The original article [6] does not claim 
this, but rather computed the co-dimension of the space of Hermitian 
matrices with degenerate eigenvalues in the whole space of Hermitian 
matrices. This is where the pure dimension count takes place. 

In our setting the calculation proceeds as follows. The function P 
on IR 3 x R is given by P(a, b, c, z) = z 2 — a 2 — b 2 — c 2 . The singularity is 
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A\ which has a one dimensional base A = R. The characteristic map is 
the map H : (a, b, c) H- —a 2 — b 2 — c 2 . The discriminant is just the point 
G R and we see that we have a level crossing over H _1 (0) = (0, 0, 0). 
And indeed, the inverse image is zero dimensional and the codimension 
of its locus is 3. All other fibers of S are of codimension 1 as one would 
expect from just a count of equations. Here we neatly see how the 
naive equation count fails over the special fiber. 

The extra dimension drop can be explained using Picard-Lefschetz 
theory as a vanishing sphere. Likewise the "diabolical" nature of these 
points — that is the behavior of wave functions when moved around 
the conical singularity — can be explained via the classical monodromy 
operator El. 

The fact that the singularity is conical can readily be checked in our 
framework. Indeed P has an isolated critical point at (0, 0, 0, 0) with 
value and signature ( h) of the Hessian. 

Finally, if one looks at the even larger family, H(a, b, c, d) = aa x + 
bo y + co z + did of all Hermitian 2x2 matrices, then one gets a family 
over 1R 4 with P(a,b,c,d,z) = z 2 - 2dz + d 2 - a 2 - b 2 - c 2 . So in 
this situation one has to use the shift s : z — > z + d upon which 
P(a, b, c, d, z) = z 2 — a 2 — b 2 — c 2 . The characteristic map 5 (a, b, c, d) = 
—a 2 — b 2 — c 2 , and we see that singularities are over S _1 (0) = (0, 0, 0, d) 
so that the fiber is now one dimensional, but still of codimension 3. The 
singular locus in the spectrum is given by a conical singularity crossed a 
line and hence not isolated. Indeed the Jacobian of P vanishes along the 
line (0, 0, 0, d, d) with critical value and the Hessian of P is degenerate 
at these points. 



2. Graph Hamiltonian and Wire network setup 

We now discuss the families that come about by considering Hamil- 
tonians obtained from finite graphs with (commutative) C* algebra 
weights on the edges. These in turn arise from wire networks of real 
materials. Here the C*-algebra in question is the noncommutative 
n-torus Tq where is a skew symmetric matrix that encodes the com- 
mutation relations of the n unitary generators Ui of Tq. Physically it 
corresponds to a constant magnetic field B. 

The initial setup for graph Hamiltonians works with a general C*- 
algebra s$ ', the relevant Hilbert space being an C 2 space. In case srf is 
commutative, by the Gel'fand-Naimark theorem, we get another de- 
scription of the algebra that srf and the Hamiltonian generate as the set 
of levels of a family of finite dimensional Hamiltonians parameterized 
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over a base. For concreteness, we will set srf = Tq, but also comment 
on how to treat the general case. 

2.1. Hamiltonian from a finite graph. In order to set up the gen- 
eral theory for graph Hamiltonians, we fix a finite graph f, a rooted 
spanning tree r of f , an order < of the vertices of T such that the root 
of r is the first vertex, a skew symmetric matrix O, and a morphism 
w : {Directed edges of T} — > Tq which satisfies the following 

(1) w(e) = w(e)* if e and e are the two orientations of an edge e. 

(2) w(e)w(e) = 1 

(3) w(e) — 1 G T e if the underlying edge e is in the spanning tree. 

Let k be the number of vertices of T. We will enumerate the vertices 
v o, • • • , v fc_i according to their order; Vq being the root. Given this data, 
the Hamiltonian H = H(T, r, <, w) G M fc (T e ) is the k x k matrix whose 
entries in T2, are 



(4) H l3 = w ^ 

directed edges e from Vi to Vj 



2.2. Family in the commutative case. If © = we can consider 
these Hamiltonians as a family of Hamiltonians over the n-torus T n as 
follows. First Tg is the C* algebra of continuous C-valued functions on 
T n via the Gel'fand-Naimark correspondence. Considered as a space 
each point of Tq is given by a character or C*-algebra morphisms 
X '■ Tq — >■ C. To each such a point, we associate the Hamiltonian 
x{H) G Mfc(C) where x is the natural lift of x to the matrix ring. 
Specifically, 

(5) (x(H)h = xiHij) 

The correspondence between points and characters is in the follow- 
ing way. Each point t G T n gives rise to the evaluation map ev(t) : 
C*(T n ) — > C which sends a function / to its value f(t) at t. Varying 
/ we obtain a character. The Gel'fand-Naimark theorem asserts that 
this is 1-1. 

Using this correspondence, we get a Hamiltonian H(t) for each t. 
Physically if we think that T n parameterizes momenta, we obtain H(t) 
by just plugging in the given momenta. 

Using the formalism we developed, we get a topological cover X — > 
T n where the points over a base point b are the eigenvalues of H{t) and 
furthermore a realization of this cover as a subspace in the manifold 
T n x R and all of our analysis applies. 
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2.2.1. C* -geometry. One can understand the topological cover X — > 
T n in C7*-geometry which yields the basic connection of our analysis 
in this section to the previous one. Consider &q C Mfc(To), the alge- 
bra generated by H G M^Tq) and the diagonal embedding of Tq into 
Mfc(Tg) as scalars. This algebra is still commutative and again by ap- 
plying the Gel'fand-Naimark theorem, we obtain a compact Hausdorff 
space X, such that is C*(X). The main point is that the cover of 
the torus given by the C* analysis from the inclusion Tq — > SS (see [I]) 
is exactly the cover 7r : X — > T n considered in the last section. 

Remark 2.1. One can readily generalize this situation to any com- 
mutative unital C* algebra srf ' . We then get a Hamiltonian over the 
base space B which satisfies C* (B) = The role of the algebra SSq is 
then played by the algebra in Mk{^) generated by H and the diagonal 
embedding of srf . 

2.3. Further characterization of P. Again considering P as a poly- 
nomial in z, the coefficient functions a k in equation ^ can be given a 
graph theoretical interpretation. For this it is convenient to introduce 
the graph T simp and the weight function w + associated to (T,w). The 
vertices of f simp are just the vertices of T. The edges of f simp are simply 
the equivalence classes of edges of T, where two edges are equivalent 
if they run between the same vertices. This identification induces a 
weight function w + on Tsimp, where now w + ([e]) = Yle'e[e] w ( e ')- That 
is the sum over all edges connecting the same two vertices as e. In this 
notation 

(6) 

1 ^ if there is no edge between w< and Vj in T smp ^ 

Hij = < w + ([e]) if there is a necessarily unique oriented edge [e] 

[ from Vi to Vj in T simp 

Plugging this into the usual determinant formula 

det(A) = ^ sign(a)a lu{1) ■ ■ ■ a ka{k) 

we can give the summand corresponding to a graph combinatorially. 
Decompose a into cycles ci, . . . , c q of length li, . . . ,l q . Then each cycle 
corresponds to a unique cycle of oriented edges in f simp . Explicitly if 
Cj = ( ji ji ■ ■ ■ jij ) then the cycle of f simp is given by the unique directed 
edges from Vj r to Vj r+1 and from Vj t , to if all these edges exist. In 

that case and if lj > 1 we set w + (cj) = T\ w + [e] where the product runs 
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over all the oriented edges in that cycle. Otherwise set w + (cj) = 0. If 
lj = 1 then w + (cj) = — z + w + $^ e (e) + w + (e) where e are loop edges 
from Vj 1 to itself if these exists and if there are no such edges, set 
w + (cj) = —z otherwise. In this notation: 

q 

(7) P(t, z) = (-l) fe sign{a)p a (t, z) with p a (t, z) = J[ w + { Cj ) 

aess k j=l 

2.3.1. Graphs with no small loops. Assume that T has no small loops, 
that is edges which return to the same vertex. Then all the diagonal 
entries Hu = and hence a&_i = 0. Furthermore, if % be the number 
of cycles of length one and we assume that these are the first i cycles, 
then 

(8) P, = (-zy n ™ + (cj)(t) 

j=i+l 

where the product is now over the cycles of length > 1. One can also 
again read off that a^-i = 0. Namely, if k — 1 cycles have length 1 then 
all k cycles have length one. 

Using the formula ((8| the coefficient a^_ 2 becomes 



(9) a fe _ 2 = - J] w + (e)w + (e) 

eer s 

/ 771 p 

2.3.2. Simply laced graphs with no small loops. Thus if furthermore T 
is simply laced, that is there is at most one edge between two vertices, 
then ak-2 = \E(T))\ is simply the number of edges. Summing up in 
this case: 

(10) P(t, z) = z k - \E f \z k - 2 + a k _ 3 (t)z k - 3 + ■■■ + a {t) 
Applying the results of §1.4.2 in this situation yields the following. 



Theorem 2.2. If the graph T has no small loops the branched cover 
X —7- T n is the pull back of the miniversal unfolding of the A^-i singu- 
larity along the characteristic map E : (ao(t), . . . , a^it))- Otherwise 
it is equivalent to the pull-back. 

The characteristic region is compact connected and lies in the closure 
of a component of the real part of W 1 ^ 1 \ D n i? n_1 over which the 
discriminant is positive. 

Moreover if T has no small loops and is simply laced, then R is 
contained in the hyperplane a k -2 = — \Er\ of A = 



SINGULARITIES, SWALLOWTAILS AND DIRAC POINTS 



17 



2.4. Wire networks. We briefly recall the relevant notions from [T] 
which we will need in the examples. As in the introduction, we fix a 
graph T embedded in M n , and a maximal translational symmetry group 
L. such that T := T/L. Let 7r : T — > F denote the projection. We will 
directly set the magnetic field form = 0. Let Vr be the set of vertices 
of T and Vp be the vertices of T and consider Jf = £ 2 {Vr) = ©„ e f ^K, 
where Jf? v = l 2 (it~ l (v)). The translation group L then naturally acts 
by translation operators on preserving the summands. This rep- 
resentation is then by commuting linearly independent unitaries and 
hence gives rise to a copy of Tq. 

Fix a spanning tree, with root vq, and an order of the vertices. Let Tg 
be the translation operator along a vector e. Notice that the oriented 
edges e lift to unique vectors in W 1 under 7r _1 . Let T ViVo = T~l be the 
total translation along the unique shortest edge path in the spanning 
tree from v to v j. Then T v . vo given an isometry J% — > Jrf? vv And hence 
we get an isometry Jt? ~ ©„ e r ^ := ^o- 

Then the weight of an oriented edge e from the vertex Vi to the ver- 
tex Vj is the translation operator T VQVi TgT ViVo . This defines the Harper 
Hamiltonian as the corresponding graph Hamiltonian which acts on 
Mq. It is shown in [I] that indeed these translation operators lie in the 
translations generated by L and hence give unitaries in Tq. Pulling 
back the translation operators of L to Mq, they together with H gen- 
erate the commutative C* algebra 

The physical background for this data as explained in detail in [TJ[2] 
is as follows. Given one of the triply-periodic CMC surfaces, P (prim- 
itive), D (diamond) or G (Gyroid), one can consider its thickened or 
"fat" version. Its boundary then consists of two non-intersecting sur- 
faces, whence the name Double Gyroid, for instance. These surfaces 
give interfaces which appear in nature. In particular, the Double Gy- 



roid could recently be synthesized on the nano-scale 12 . The structure 



contains three components, the "fat" surface or wall and two channels. 



Urade et al. 12 have also demonstrated a nanofabrication technique 
in which the channels are filled with a metal, while the silica wall can 
be either left in place or removed. This yields two wire networks, one 
in each channel. The graph we consider and call Gyroid graph is the 
skeletal graph of one of these channels. The P, D, and G examples 
are the unique such surfaces where the skeletal graph is symmetric and 
self-dual. The graph Hamiltonian is then the Harper Hamiltonian for 
one channel of this wire network. The 2d-analogue of this structure is 
the honeycomb lattice underlying graphene. 
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e 2 4 -5 5 

Figure 2. Spanning tree and characteristic region for 
the Gyroid (solid region). The curve is the slice of the 
discriminant of the A% singularity at 02 = —6 



3. Calculations 

Since all calculations are for the base T n we will consider the func- 
tion P locally pulled back via the exponential map exp : R n — > T n 
(di, . . . , a k ) — >■ (exp(zcii), . . . , exp(za„)). In this notation given a point 
(dj) and its corresponding character x, the translation operators Uj cor- 
responding to the generators of the algebra Tq get mapped to x{Uj) = 
exp(iaj). Indeed under the Gel'fand-Naimark correspondence the op- 
erator Uj is the function exp(iaj) for the coordinates of T n above. To 
simplify the calculation, we will drop the x an d just write the function 
for the operator. 



3.1. The Gyroid. 



3.1.1. The matrix and the function P. As shown in [T], the relevant 
graph for the Gyroid wire network is the full square which is simply 
laced. We also fix a spanning tree shown in Figure [2] ordering the 
vertices as indicated, the root being 1. 

The Harper Hamiltonian is a function on T 3 = S 1 x S 1 x S 1 . We 
can visualize T 3 as a cube where opposite sides are identified. 

The Harper Hamiltonian reads [l] 
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1 1 1 \ 

A B* 
A* C 
B C* J 

where A, B and C are operators generating T 3 , each of which we can 
think of as a function on S 1 We will rewrite them as A = exp(ia), 
B = exp(ib), C = exp(ic) with a, b, c real. 

The eigenvalues of H are given by the roots of the characteristic 
polynomial: 

(12) P(a, b, c, z) = z 4 — 6z 2 + ai(a, b, c)z + a (a, b, c) 
where 

ai = —2 cos(a) — 2 cos(o) — 2 cos(c) — 2 cos(a + b + c) 
a = 3 — 2 cos(a + 6) — 2 cos(6 + c) — 2 cos(a + c) 

give the characteristic map H := (ao,ai) : T 3 — )■ M 2 

3.1.2. y4 gwicA; /ooA; i/ie characteristic region. The characteristic re- 
gion of the full square graph is depicted in Figure [2] The curve shown 
in the figure is the discriminant locus which is explicitly given by 

(13) 20736 a - 4608 a 2 , + 256 a 3 + 864 a\ - 864 a a\ - 27 a\ = 

The boundaries of the characteristic region are obtained as the col- 
lection of points (ao, cli) for a = b = c and a = b = —c. 

We see that the characteristic region is contained in the slice a2 = —6 
of the A3 singularity and intersects the discriminant in exactly three 
isolated points, the two cusps and the double point of that slice of 
the swallowtail. The two cusps are in the stratum of type A 2 and the 
double point is in the stratum of type (A 1 , Ai). As is quickly seen and 
we calculate below the fibers over all these points are indeed discrete. 
For the A 2 singularities, this is just one point each, giving rise to two 
triple crossings, and the fiber over (Ax, Ax) consists of two points. Over 
each of these points there are two double crossings and it turns out, 
see below, that these are Dirac points. 

This is a very special situation in that the points on the discriminant 
are actually at singular points of the region. From the singularities it 
is easily seen that the image of the tangent spaces at the points in the 
fiber is and hence J a vanishes. 
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3.1.3. Classification of the critical points. We first check the condition 
VP = 0. This yields the equations: 

dP 

— = z(2sm(a + b + c) + 2sin(a)) + 2sin(a + b) + 2sin(a + c) = 
da 

dP 

— = z(2sin(a + 6 + c) + 2sin(6)) + 2sin(a + 6) + 2sin(6 + c) = 
ob 

dP 

— = z(2sm(a + b + c) + 2sin(c)) + 2sin(a + c) + 2sin(6 + c) = 
ac 

dP 

= _2cos(a + 6 + c) - 2cos(a) - 2cos(6) - 2cos(c) + 4z 3 - 12 z = 

oz 

To solve these equations, we rewrite the first three using trigonomet- 
ric identities as 

{ b + c r „ 

2COS( — - — ) = — cos( — - — ) 

z cos( ^ ) = — cos( — - — ) 

,a + b 
zcos{ — - — ) = — cos( 

In the case that all cosines are different from zero, we can solve each 
equation for z and set them equal. This leads to 

b — a + 2c a — b + 2c 
cos( ) = cos( ) 

2a — b + c 2a + b — c 
cos( ) = cos( ) 

From the first of these equations we get a = b mod 27r, from the 
second b = c mod 2n, and zcos(a) = —1 for cos(a) ^ 0. Plugging this 
back into the last of the original equations, we find 

8cos 6 (a) + 4- 12cos 2 (a) = 

Among the solutions we pick those for which the characteristic poly- 
nomial P(a,b,c,z) (see Eq. (12)) is zero, namely cos(a) = —1/z = ±1, 
i.e. a = 0, 7r and z = ±1. 

In the case that cos(a) = we obtain 

4z 3 - 12z = 

which has solutions z — 0, ±v3. z = has to be discarded since it 
does not satisfy P(a, b, c, z) = 0. 
Summing up, the critical points are 

(1) a = b = c = (mod 2tt); z = —1 

(2) a = b = c = 7r (mod 2ti); z = 1 
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(3) a 



2' 2 



(mod 27r); z = ±\/3 



Looking at the image of these points under the characteristic map, we 
see that 2(0, 0, 0) = (—9, —3) is the lower cusp, which is an A 2 point 
with a triple degeneracy, so is a(7i,7i,n) = (9,3). For the other two 
points, a(ir/2, tt/2, tt/2) = a(37r/2, 3tt/2, 3tt/2) = (0, 9) and they are in 
the (Ai,A\) stratum. This means that these points are candidates for 
Dirac points, which they indeed are. 

To decide this, we calculate the Hessian. 



Plugging in a = b = c 



3- 



(14) 



Hess 



2 ' 

/ 



V 



mod 27r; z = ±a/3, it becomes 
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-2 
-2 
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/ 



which has signature( h). Notice that the corresponding cone 

is also not tilted. For the other two points, the Hessian vanishes, as 
expected. 

Since the singular points are all isolated, we can take any direction as 
transversal. It turns out that the diagonal curve C given by a = b = c 
gives a transversal direction for all 4 singular points at once. The 
spectrum on this line is given in Figure [3] can be obtained by using 
group theory cf. |8j. This is how it was first found in [7] where the 
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A 

Figure 4. The graphs f for the diamond (left) and the 
honeycomb case (right) 

authors considered an equivalent family in another context. They also 
found no other singularities numerically. Explicitly the spectrum along 
this diagonal is given by 

Ai = co expiia) + Co exp(-ia) \ 2 = u exp(ia) + uj exp(za) 
(15) A 3)4 = cos(a) ± a/ cos 2 (a) + 3 

From this one can see a linear dispersion relation in the direction of the 
diagonal. Without further analysis, one cannot deduce the dispersion 
relation in any other direction from this results. 

By our previous analysis we have, however, proven analytically, that 
there are indeed no other singularities and furthermore have deter- 
mined that there is a linear dispersion relation in all directions at the 
points (tt/2, 7r/2, tt/2) and (37r/2, 37r/2, 37r/2) establishing that these 
are indeed two Dirac points. 

It is interesting to note that the image curve H(C) runs first from 
(—9, —3) to (0, 9) along the boundary of the region, continues as the 
boundary curve to (9, —3), and then turns back on itself to cover these 
two boundary pieces twice. 

The nature of the two triple points is that they are isolated and they 
are a pullback of the unfolding of the A 2 singularity. This is given by 
the image of the characteristic map where now we consider the local 
neighborhood of the cusp point in the slice as a miniversal unfolding 
of A 2 . This is indeed possible and a standard way of embedding the 
unfolding of A 2 into that of A3 [4]. 

3.2. The honeycomb and diamond cases. We treat the diamond 
and the honeycomb case in parallel. The graphs T are given in Figure 

m 
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The Hamiltonians are 

1+U+V 
1 + U* + V* 



(16) Hhon 

and 

(17) H D -- 



1 + U + V + W 

1 + U* + V* + w* 

We again use U = exp(iu), V = exp(iv), W = exp(iw). 

The polynomials are P(u, v, z) = z 2 — 3 — 2cos(u) — 2cos(v)— 2cos(m — 
v) and P(u, v,w,z) = z 2 — 4 — 2 cos(u) — 2 cos(f ) — 2 cos(w) — 2 cos(u — 
v) — 2cos(m — w) — 2cos(f — u>). The characteristic regions in R are 
just the intervals [—9,0] and [—16,0]. The discriminant is the point 0. 
From this we see that in both cases we have to have ao = and the 
singular locus is simply this fiber. 

3.2.1. The honeycomb case. For the honeycomb, the standard calcula- 
tion shows that in this case U = V* and U G {p^ : exp(27rz/3), p 3 }, 
which means that the fiber consists of 2 points. These are the well 
known Dirac points (ps, pz), (p~s, ps). We can check this explicitly: 
V(P) = (2sin(w) + 2sin(w — f),2sin(f) — 2sin(w — v),z) from which 
we see that z = and u = — v = 2-u(27r). Furthermore 

(18) 

/2 cos(-u) + 2cos(u — v) —2cos(u — v) N 
Hesshon = — 2(cos(u — v )) 2 cos(v ) — 2 cos(m — v ) 
V 1, 

which has the correct signature ( h) at the given points, from which 

we recover the known result that these points are Dirac points. 

3.2.2. The diamond case. The equation for the fiber over 0, 

— 4— 2cos(u)— 2 cos(w)— 2 cos(w)— 2 cos(w— v )— 2 cos(u— w)— 2 cos(v—w) ■ 

has been solved in [2] and the solutions are given by (u, v, w) = ((pi, <pj,4>k 
with (pi = 7T, <j>j = <pk + 7T mod 27r with {i, j, fc} = {1, 2, 3}. So in this 
case the fiber of the characteristic map is 1-dimensional and the pull- 
back has singularities along a locus of dimension 1, which also implies 
that there are no Dirac points. Geometrically the singular locus are 
three circles pairwise intersecting in a point. 

3.2.3. The characteristic map and region. For both the honeycomb and 
the diamond graph, the relevant singularity is A\. Both these graphs 
are not simply laced, so their image is not contained in a slice. The 
swallowtail is only one point and this is the stratum of type A\. In the 
honeycomb case the fiber over this point is discrete and consists of two 
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Figure 5. Spanning tree and characteristic line for a 
triangle with single bonds 

points, while in the case of the diamond lattice the fiber is not discrete 
and it is given by three circles pairwise intersecting at a point. It turns 
out that in the honeycomb case the two candidates for Dirac points are 
indeed Dirac points. While for the D case there is a non-trivial fiber 
which is essential 1-dimensional. Hence we do not get Dirac points, 
but rather spread out singularities. 

3.3. Three— vertex graphs. In order to have some examples for the 
A2 singularity and to show the kind of behavior that is possible, we 
considered a three vertex graph with either only simple edges, or one, 
two or all of the edges doubled. 

The characteristic regions are seen in Figures [5](9j Here we see that 
for the simply laced case, we get a slice, for one or two doubled edges, we 
get parabolic regions, which intersect the boundary in two points, which 
are of type A\ and finally, in the case of all edges being double, the 
characteristic region is a surface, which is bounded by the discriminant 
and the line on which a\ takes its maximal value, in this case a± = 12. 

3.3.1. Triangle with single bonds. We consider the graph and the span- 
ning tree given in Figure |5j 

The associated Harper Hamiltonian reads: 



where A is an operators on S 1 . We will rewrite it as A = exp(m) with 
a real. The characteristic polynomial is: 



The characteristic region is easy to calculate, since the graph is 
simply laced it is contained in the slice a\ = —3. The image under 



(19) 




(20) 
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-2cos(a) = [-2,2] thus R = [-2,2] x -3 in IR 2 . Figure § shows this 
together with the zero locus of the discriminant. 

From this we see that all possible singularities occur at ao = —2 cos(a) 
±2 that is a = 0,7r(27r). Indeed calculating VP(a,z), we use 

(21) f a = 2sin(a) 

(22) - = 3*»-3 

These equations vanish simultaneously for the following choice of 
variables: 

(1) z = ±l 

(2) a = 0, 7r (mod 2tt) 

Among all possible combinations, the choices (z = —l,a = 0) and 
(z = 1, a = 7r) are zeros of p(z, a). For these two points, the Hessian 

(23) He 8S =( 2c< f> I 

has a no n- vanishing determinant of det(Hess) = —12 and signature 
( — h). So we find two Dirac points. Again the cone is not tilted. 



3.3.2. Triangle with one double bond. We consider the graph and the 
spanning tree given in Figure [6] where one of the bonds is a double 
bond. The Harper Hamiltonian reads in this case 

/ 1 1 

(24) H= 1 A + B 

\ 1 A* + B* 

where A, B are operators on S 1 . We will rewrite them as A = exp(ia), 
B = exp(ib) with a, b real. The characteristic polynomial is: 

(25) P{a, b, z) = z 3 - (4 + 2 cos(a - b))z - 2 cos(a) - 2 cos(6) 

Again, we set a\ = —(4 + 2cos(a — b)) and ao = — 2cos(a) — 2cos(6). 
We see that this time the characteristic region is not contained in a 
slice, which was not to be expected since the graph is not simply laced. 
The region is depicted via a scatter plot in Figure [6} One reads off that 
R intersects with the discriminant locus in two points. 
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Figure 6. Spanning tree and characteristic region with 
double bond 



To calculate VP(a, b, z) = we use, 

dP 

da 
dP 

~db 
dP 

dz 

From the first two equations, we get either a 



(26) 



2 sin(a — b)z + 2 sin(a) 
—2 sin(a — b)z + 2 sin(6) 
(4 + 2cos(a-6)) 



3- 



0, 7r and 6 = 0, 7r, but 



for all combinations of those, the remaining two equations (ff = 
and P(a, b, z) = 0) cannot be simultaneously solved. Therefore the 
only possible solution for the first two equations is to take a = —b and 
z = — sin(2a) f° r s i n (2a) 7^ 0. Putting this into the last equation yields 
the trigonometric equation 

3 sin 2 (a) - 4 sin 2 (4a) - 2 cos(2a) sin 2 (2a) = 

which has the solutions a = ±|, ±^p These also lead to a vanishing 
of the characteristic polynomial. 

So we get the following two solutions (the negative values lead to the 
same values for ao and ai): 



3 ' 
2tt 



-| mod 27r, z = — 1 
mod 2tt,^ = 1 



(1) a 

(2) a- ,.. 

The Hessian is 
(27) 

2 cos(a — b)z + 2 cos(a) —2zcos(a — b) 2sin(a — b) 

Hess = I — 2zcos(a — b) 2 cos(a — b)z + 2 cos(6) — 2sin(a — 6) 

2sin(a — 6) — 2sin(a — 6) 6z 

and has signatures (+ H — ) and ( h), respectively. Here the cone 

is actually tilted. Again there are two Dirac points. 
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3.3.3. More variations. In the same way, we can obtain information 
about possible Dirac points for the following graphs: 

al 



-6 -4 -2 J 


\ 2 4 6 


A 

MM 

-12 


■At 

h 



Figure 7. Spanning tree and scatter plot of character- 
istic region with triple bond 

al 




Figure 8. Spanning tree and characteristic region with 
two double bonds 



al 




Figure 9. Spanning tree and scatter plot of character- 
istic region with three double bonds 



For the triple bond case shown in Figure [7] and the two double bonds 
in Figure [8j we see two isolated intersection points, but the fiber will 
be dimension 1, so like in the D-case there will be no Dirac points. 
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Considering three double bonds (see Figure [9]) we see that now the 
intersection of R with D is along two of the boundaries of R and hence 
not only will the fiber have dimension 2, but we also expect to have 
horizontal directions which to not resolve the singularity. We leave this 
for further investigation. 

3.4. The P and other Bravais cases. Here the graph has small 
loops, and the characteristic polynomial has to be transformed. It is 
simply a polynomial of degree 1. P(t, z) = z — J2i U, so that after shift- 
ing z we are left with just z — 0, which is not critical. Not surprisingly, 
there are no singularities. 

4. Conclusion and outlook 

We have developed a general method to analyze singularities in the 
spectra of smooth families of k x k Hamiltonians parameterized by a 
base B using singularity theory. In particular, we realized the spectrum 
as a singular submanifold X of the smooth product BxRas the zero set 
of a function P. This led us to a simple characterization of Dirac points 
as A\ or Morse singularities of P with critical value whose signature is 

(— •• h) or (+ • • • H — ) . Furthermore we could represent tt : X — > B 

up to a given diffeomorphism on the ambient space as the pull-back 
of the miniversal unfolding of the A^-i singularity via a characteristic 
map S. This classifies all the possible singularities of the fibers of 
it as (An, • • • j AiJ The image of the characteristic map, called the 
characteristic region, allows one to read off which ones occur. 

We then applied these techniques to the graph Hamiltonians and wire 
networks. Here we reproduce the known Dirac points for graphene and 
make the surprising find that the Gyroid wire network also has Dirac 
points. We expect that this should have practical applications. 

The situation for the Gyroid is very special, as the characteristic 
region goes into the cusps and the self-intersection locus of the swal- 
lowtail without prior contact to the "walls" . Were this not be the case, 
one would not expect isolated points. We gave more graph examples 
to illustrate how special this behavior is. Adding multiple edges, we 
expect to get "full" region as in case of the three double bonds in a 
triangle. 

One exciting find is that there seems to be a commutative/non- 
commutative duality in the wire network families first stated in [2] . By 
this we mean the observation that there is a correspondence between 
the locus of degenerate points in the commutative setting and the locus 
of parameters where the corresponding non-commutative algebra 38® 
is not the full matrix algebra Mfc(Te). The correspondence is not 1-1, 
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but the top dimensions agree and there are further features that look 
dual. We wish to emphasize that although the space T n appears in both 
settings as the parameter space, it a priori plays two different roles. 
In the case without magnetic field, the parameters are quasimomenta, 
while in the case with magnetic field they are the field's components. 
It is intriguing to speculate that the non-commutative setting is a 
model for a non-commutative unfolding of singularities and that this 
furnishes the framework to make the duality explicit. This will be a 
topic of further research. 

There are several other directions of research that are immediate. 
First one can ask if in the graph case there are symmetries forcing 
the degeneracies. This is pursued in [8] using a regauging groupoid 
action. Second, a physically relevant question is how stable the singular 
fibers are with respect to deformations of the Hamiltonian. For a 3- 
dimensional parameter space, simple singularities, i.e., Dirac points, 
are "magnetic monopoles" [3] and are expected to be topologically 
stable. This, and the evolution of other types of singularities, is further 
discussed in |13|. We will also focus on making the type of analysis 
explicit in the non-commutative geometry language. There should be 
some kind of characteristic classes and parings much like in the setting 



of the quantum Hall effect as presented in 14, 15 
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